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ity field equations. Outside the (delta-like) source this not a vacuum so- 
^jrj' lution of General Relativity. It represents the contribution of the massive, 

m = (— /3)~^/^, spin-two field associated to the Ricci squared term in the 
gravitational Lagrangian. The fourth order terms tend to make milder the 
singularity of the curvature at the point where the particle is located. We 
\ generalize this analysis to L'-dimensions, extended sources, and higher than 

fourth order theories. We also briefly discuss the scattering of fields by this 
kind of plane gravitational waves. 
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I. INTRODUCTION 



Higher order theories of gravity are the generally covariant extension of General Rela- 
tivity when one considers in the Lagrangian nonlinear terms in the curvature. The field 
equations derived by variation of this Lagrangian contain derivatives of the metric of an 
order higher than the second (except for the Lovelock's Lagrangian constructed from the 
D-dimensional extension of four dimensional invariants [|l|). Historically, they have been 
introduced by Weyl right after General Relativity |^ . These theories provided a framework 
where to study the unification of gravity with other fundamental fields, and the possibil- 
ity of classically avoiding cosmological singularities. More recently, higher order theories 
have been shown to lead to inflation and dimensional reduction without the introduction 
of any additional scalar field. Nowadays, among the main motivations for their study are 
their appearance (as vacuum polarization terms) in the one-loop renormalization of fields in 
curved spacetimes and in the low-energy limit of string theory ||^. 

In four dimensions, using the Gauss-Bonnet invariant, a general fourth-order (quadratic) 
covariant Lagrangian can be written as 



I d^x I - 2A + + ^ PRf.^R^"' + 167rG£„| , (1.1) 



where we have not considered surface terms since they will not contribute to the analysis of 
the field equations we will perform. 

The field equations derived by extremizing the action / are given by 

X 1 G vi Cj dl 

Rfiu - -^Rgfiu + J^Qfiu + oiH^y + /?J^j, = STrGT^i, = -—j==-^^ , (1.2) 



where 



and 



H^u — —2R.^y + 2g^yUR — -g^^R'^ + 2RR^i, , (1.3) 



1 _„ 1 



Ifiu — — 2_R|J .^^ + D-R^j. + -g^^UR + 2R°Rau — -g/iuRajsR"''^ ■ (1-4) 



Upon linearization of these equations one can see that in addition to the usual graviton 
field this theory contains a massive scalar field related to R and a massive spin-two field 
ipfj^u related to Rf^i, (see Ref. for further details). Asking for both fields to have a real 
mass (to recover the Newtonian limit) leads to the, so called, no-tachyon constraints 

3a + p>0 , p<0. (1.5) 

The value of these coupling constants can only be determined by experiments or could be 
computed from a fundamental theory that would unify gravity to the other forces in nature. 
It is then expected they to be of the order of the Planck scale. 

The quantum properties of these theories have extensively been discussed However, 
several issues, such as the unitarity problem p[ and the semiclassical instabilities ||^, remain 
to be solved. Besides, we have not even a complete understanding of the proper solutions 



2 



of fourth-order theories. In Refs. |T0| we have developed a perturbative method to find 
solutions of the field equations (|1.2|) , given a solution to the general relativistic problem (see 
also Refs. |TT|,p|). The method essentially consists of writing the higher derivative terms as 
derivatives of the matter energy- momentum tensor T^j,. One obtains a series development 
around the General Relativity metric in powers of the coupling constants a and (5. We have 
thus studied the a and /? corrections to the Reissner-Nordstrom and straight cosmic string 
metrics P^ . 

Our perturbative approach, evidently, gives no corrections when the matter energy- 
momentum tensor vanishes. In other words, the method confirms that vacuum solutions 
(including A 7^ 0) of General Relativity are also solutions of higher order theories. However, 
the converse is in general not true. Higher order theories have a richer set of vacuum 
solutions than General relativity. If we call this set T^vho and the set of vacuum solutions 
of General relativity T,vgr, its difference, Sav = ^vho ~ ^vgr, is in general a non-empty 
set. From the above described inability of the perturbative approach to find solutions in 
Eav, we can deduce that such solutions, if they exist, have to be non-perturbative around 
a and /5 equal to zero. 

For black hole [0 and de Sitter |[T^ solutions one can extend the no-hair theorems valid 
for General Relativity to fourth order theories in the case (3 = 0. From where we can infer 
that Sav black holes will have the form of the Kerr metric plus non-analytic corrections in 
(3 only. The perturbative corrections in powers of (3 to the Reissner-Nordstrom metric have 
been given in Refs. |T5| , p!0[| . 

In Sec. II we introduce the gravitational shock waves, a special case of pp-waves solution 
of Einstein equations with a delta like source term. We analyze its extension to fourth 
order gravity and find the corresponding solution which is an exact proper solution of field 
equations ( |1.2D . In Sec. HI we generalize this exact solution to D dimensional spacetime, 
extended sources, and theories of order higher than the fourth. We finally, in Appendix A, 
deal with the problem of the scattering of a scalar field by these shock wave geometries and 
compute the S'-matrix for the case of a source obtained by boosting to the speed of light 
the Kerr metric. 



II. GRAVITATIONAL SHOCK WAVES 

As pointed out by 't Hooft [|16|, at energies of the order or higher than the Planck scale 
the picture of particles propagating in flat spacetime ceases to be a good approximation. 
Interestingly enough the curved metric generated by such particles has a remarkable simple 
form 

ds'^ = —du dv + H{u, x±) dv? + dx]_ , (2.1) 

where 

H{u,Xi_) = f\xA_) 5{u) , u = t-z, v = t + z. (2.2) 

This metric represents an impulsive gravitational wave localized in the plane u = 0, i.e. along 
the motion of the particle. The shock wave is accompanying the particle, both traveling 
at the speed of light. The profile function f{x±) is the only quantity depending on the 
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characteristic of the source. It is only a function of the coordinates along to the plane of 
motion, x^. The geometry is just flat before and after the pass of the wave, i.e. for m 7^ 0, 
and is an special case of plane fronted with parallel rays (pp) - waves |T^. Geodesies are 



then just straight lines outside the wave front and have a finite discontinuity (or shift) in the 



V coordinate as they cross u = given by (see Ref. |I8| for the case of the Aichelburg-Sexl 
metric, i.e. Eq. ( |2.10|) below) 

At; = -f{pi) , (2.3) 

where pi is the coordinate distance from the origin (where the source is located) to the point 
where the geodesic reaches m = 0. 

There is also an effect of spatial refraction of geodesies (see again Ref. 

1 

— ( 
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cot(^,0 + cot(ft,e) = 7;dpf{Pi) , (2.4) 



where 6'j„ and Ore are the incident and refracted angles respectively. Of course, locally 
physical measurable quantities involve the relative shift or refraction of nearby geodesies 
which involve one further derivative with respect to p, the cylindrical coordinate on the 
plane u = 0. We shall came back to this point by the end of this section when we will 
compute the Riemann tensor components which give the geodesic deviation. 

The only non-vanishing components of the Riemann tensor for metric ( |2.1j ) are (apart 
from the ones obtained by symmetry properties) 

Riuju = -^dijH{u,x_L) , i,j = . (2.5) 

And the only non-vanishing components of the Ricci tensor for metric ( ^2.11 ) are (apart 
from the ones obtained by symmetry properties) 

Ruu = ~VlH{u,x^) , (2.6) 

where V^^ stands for the Laplacian operator in the Xj_ space. 

The curvature scalars formed out of the Ricci squared, Kretschmann, and Curvature 
scalar tensors all vanish identically since, 

R = 0, Raf,R'"' = 0, Rap^sR''^^' = . (2.7) 

To determine the form of the profile function one considers a null source represented by 
the following energy-momentum tensor 

Tun = (t{x^)6{u) , (2.8) 

with all other components vanishing. 

Next, one has to impose the field equations to the metric. Einstein equations linearize 



and reduce to a Poisson equation in the x± space. [lT9| , p!8 

1 



^S{u)Vlf{x^) = MGT^u = 87cGa{x^)6{u) . (2.9) 
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The simplest source one can consider is a spinless particle, with momentum p, represented 
by ap{x±) = pS{x±) . Then, the solution of Eq. ( p.9|) is readily found to be (in four dimensions 
for the sake of simplicity) 



(2.10) 



where now the profile function only depends on the radial distance to the origin, p, in 
the plane u = (as expected for a cylindrically symmetric problem), po is an integration 
constant having the units of p. 

Metric ( |2.1| ) and ( |2.10| ) was originally obtained by Aichelburg and Sexl by boosting 
Schwarzschild metric along the z-axis and taking simultaneously the limits M ^ and 

The explicit form of the profile function have also been found for a variety of sources 
2l|^5| by use of the above sketched procedure. 

As explained above one expects the shock wave metric to be relevant in processes in- 
volving energies of the order or higher than the Planck scale. At such huge energies, higher 
order corrections to the gravitational theory will also be relevant. It is of interest, then, to 
study the form of the profile function, /, in the fourth order theory of gravity ( |1 . 1|) . 

If one plugs the ansatz ( |2.1| ) for the shock wave into the fourth order field equations 
( |1.2|) , obtains again a linearized equation for the profile function 



/5Vi + vl H{u, x^) = MGTu 



--5{u) [/?Vi + Vi] /(xx) = 87rG(T(xx)(5(M) . 



(2.11) 



where V^ = V^V^ is the fourth order Laplacian in the perpendicular space. 
This equation can be integrated twice to give 

/3Vl + ll H{u,xa_) = Hgr{u,x^) . 



(2.12) 



We have thus reduced the original fourth order problem to a second order one. Now, by 
making the decomposition 



H{u, x±) = Hq{u, x_l) + Hgr{u, x±) 



(2.13) 



where the index Q refers to the purely quadratic part of the solution, and using the fact 
that Hgr satisfies general relativity equations, Eq. ( p^.l2D can be re-written as 



V + 



Hq(u,x_l) = IGirGTu 



(2.14) 



We first note here that the above form of equation (but with different constant coefficients) 
have been found when analyzing shock waves on curved backgrounds in pure General Rela- 
tivity [|IB],^BH2§. 

Let us now consider the simple example of the particle represented by the source term 
o'p{x±) = p6{x±), (in D = 4). The left hand side of Eq. (|2.14| ) takes the form of a Bessel 
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equation. Our problem can thus be solved in terms of the Green function for a cylindrically 
symmetric problem. The solution can then be expressed as 



Hq{u,x^) = fQ{p)6{u) , fgip) = -8GpIo{0)Ko 



(2.15) 



where Iq and Kq are the modified Bessel functions and where we have taken into account 



that (3 < from the non-tachyon constraints (|1.5| ). 

Finally, the profile function generated by a massless uncharged particle is 



/(p) = -8Gp 



P 



+ In 




(2.16) 



Equations ( |2.1|) and ( p.l6| ) thus represent the exact solution to our problem in fourth 
order gravity. Note that the coupling constant a does not appear in the solution. This can 
be traced back to the fact that the scalar curvature R, identically vanishes for metric ( p.lD , 
as can be readily verified. 

Since in the asymptotic regime, for both small and large values of p, 



In 



2p 



-1/2 



, for small p/ y/—l3 

exp ^ T^^ ^ , for large p/ \J —fi 



(2.17) 



we have then computed the contribution to the shock wave of the massive (with mass 
m = 1/a/— /?) spin-two field. In the case of small p, the profile function assumes the form 
of an Aichelburg-Sexl profile for a point-like source. Also note that the dependence on /5 is 
clearly non-perturbative, as expected from our comments in the introduction. In fact, for 
p 7^ 0, the profile ( |2.16| ) gives a vacuum solution to the fourth order theory which is not 
solution of General Relativity, i.e. it is in the set Sav- For /? — we recover the general 
relativistic metric; while p = (no source) gives the flat space limitfj. 

In Fig. |1] we plot the resulting profile function and compare it to the general relativistic 
one. The first apparent difference is the non divergence in the origin of coordinates, i.e. where 
is located the particle. Even if, as we have seen in (Eq. ( p.7| ), for plane gravitational waves 
the curvature scalars identically vanish, this does not necessarily imply the gravitational 
field will be non-singular at the origin. In fact, for the cylindrically symmetric situation we 
are analyzing the only non-vanishing components of the Riemann tensor for metric ( p.l|) 
are (apart from the ones obtained by symmetry properties) 



R 



-^dlH{u,p) , R^ucpu = -^dpH{u,p) 



pupu 



(2.18) 



^Note that here we refer to the no source case as imposing a boundary condition on the solution 
such that it represents Minkowski space; while, when we refer to a vacuum solution for p 7^ 0, it is 
the result of imposing other boundary conditions such that Eq. ( 2.16| ) is the solution. It is in the 
same sense that one refers to the Schwarzschild solution as a vacuum solution of Einstein theory, 
while for M ^ one finds Minkowski space. 
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and the only non-vanishing components of the Ricci tensor for metric (|2.1j ) is 



Gun = ~Yp^p ip9pH{u,p)) = ^Hq{u,p) 



It then follows that components of the Riemann tensor like 

1 „ . 1 



Rl,=d'H{u,p), Rl 



dlH{u,p) 



(2.19) 



(2.20) 



will diverge logarithmically as p — 0. This can be classified as a "parallelly propagated 
curvature singularity". However, we note here that the singularity is notably milder than 
the 1/p^ for General Relativity. 



III. GENERALIZATIONS 



Many of the present candidates to unify gravity with other interactions consider D, the 
dimensionality of the spacetime, bigger than four. Since Eq. (|2.11|) can be extended to any 
dimension, we find that the D-dimensional generalization of Eq. (|2.1(j|) is given by 



/(P) 



levrGp 



D-3 



-2f3) 



2-D/2 



2~D/2 



A~D' 



T{D/2-l] 



K. 



2-D/2 



-(3 J {4-D)\po, 



(3.1) 



where Qd-s = 2n^/^-^ /T{D/2 - 1) is the area unit in the D — 3 sphere and Gp carry units 
of length^''^ . 

In the case the source term is extended, but keeps its axial symmetry, i.e. a{xi_) = a{p) 
we have 

/(p) = fGR{p) + 167rGx2-^/2 r {K2-D/2{x)lD/2-2{r) " I D /2-2{x)K2-D /2{r)] r""^'^ a{r)dr 

(3.2) 

where as before the index GR reffers to the solution of the problem in Einstein theory and 
X = p/ \/ —(5. When the source has a Q dependence, we can Fourier transform it as well as 
the solution and it will look like Eq. ( |3.2| ) with the index of the Bessel functions being I now 
(in D = 4 and where / refers to the corresponding Fourier mode in Q). 

An interesting example of extended source easily solvable is the boosted straight string 
2^ . In four dimensions Eq. ( |2.11| ) becomes an ordinary fourth order differential equation 



with constant coefficients and a 5(y)-like source. The profile function takes then the form 



f{y) = -8nGp 



-(3e 



+ \y\ 



(3.3) 



where y refers to the perpendicular distance to the string measured on the plane ti = 0. One 
can check that this exactly corresponds to a particle in D = 3 dimensions from expression 

One can also see the problem of quadratic theories as giving a correction to the energy- 
momentum tensor (as in the semiclassical approach to quantum gravity) and think of the 
problem as one for Einstein gravity with an effective source. Then, from Eq. (|2.12|) we have 
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^Vi/(x^) = SnGaeffix^) = ^(/ - Jgr) ■ (3.4) 



Hence the explicit form of (Te// in terms of integrals of cr can be read off of Eq. (|3.2|) . For 
example, for the point-like source considered in the generalization of the Aichelburg-Sexl 
metric, we obtain from Eq. ( p.l|) 



^]z)-3 r(D/2 - 1 




-e//(p) = TT^ ^Tr. ( ^ 1 K,-o,. ( ^ 1 , (3.5) 



which represents an extended effective source. 

Another kind of generalization is to consider a gravitational Lagrangian containing terms 
in the curvature and its derivatives that will generate field equations with derivatives higher 
than the fourth, for example, terms like i?^, RR^^R'^^, R^y-\R^''''^ , etc. When we consider 
solutions to the field equations of the form Eq. ( p.l|) there will not be any contribution 
coming from the terms i?" with n > 2 since i? = for this metric. Neither terms with 
contractions of the curvature tensors involving powers higher than the second, since the only 
non vanishing component of the Ricci tensor is Ruu- Term involving covariant derivatives 
of the curvature tensors (like R^^^xR^'^''^) will, however, give a contribution in the form of 
higher order Laplacian operators (like Vj^ in the above example). For a theory containing 
terms of up to power N in the curvature tensor and its derivatives we expect to have the 
following form of the field equations for a metric of the type ( [^.ID 



1 ^ 1 



Q Z I / " ■ _L \~i ~ r\ 



N 



n V 



.n=l 



72 2 



H{u,x±) = SttGTuu , (3.6) 



where the m„'s are constants related algebraically to the 7„'s and can be seen as the masses 
of some of the particles of the theory, for example, mi = corresponds to the graviton, 

= to the massive spin 2 field of quadratic theories, and so on. 

Let Hn be a solution of the equation (V^ —m'^)Hn = 0. Then Hhom = Y.n=i CnHn (with 
Cn constants to be determined by imposing the boundary conditions of our problem), will 
be a general solution to the homogeneous problem associated to Eq. ( |3.(j| ) as can be easily 
verified. We thus construct the general solution to the inhomogeneous problem by adding a 
particular solution that can be found by use of the Wronskians method, i.e. 

HM^p) = j:H^{u,p)f'^^dr. (3.7) 



n=l 



Still a further generalization of the solution can be done by considering the wave traveling 
on a curved background instead of a flat one. This can be performed following the steps Dray 



and 't Hooft made in Ref. [T3] for Einstein theory and should present no further difficulties 
than obtaining the coefficients in the resulting field equation at m = (generalization of Eq. 



IV. DISCUSSION 



We have seen that a nice feature of the exact solution Eq. ( p.l6|) is the fact that it makes 



milder the curvature singularity at the location of the particle generating the gravitational 
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field. It is in fact one of the historical classical motivations to introduce higher order theories 
in Cosmological scenarios as we recalled in the introduction. It is also a desirable feature 
when one thinks of this quadratic theory as an effective one being the by-product of the 
low-energy limit of a finite quantum theory of gravity. We have studied an exact solution 
of the theory given by the action (|1 . 1| ) . Among the motivations for studying this kind of 
theory we gave the argument of its similarity with what one finds renormalizing to one-loop 
a field theory in curved backgrounds. In this case one has also to take into account non-local 
corrections in the renormalized energy-momentum tensor. That, in fact, can be done taking 
into account the results of paper Let us recall that the final result in this case would be 



only valid to order h since the renormalized energy-momentum tensor was computed to that 
order. This point bring us to the question of "self-consistency" of the perturbative approach 
discussed in Ref. 0. There it is in general considered only solutions linear in h (and in the 
coupling constants a and /3) for consistency with the field equations which are considered to 
be only precise to one-loop. This point of view has the advantage of avoiding the, so called, 
runaway solutions; rendering thus, Minkowski space stable. In our view this procedure is too 
restrictive and precludes some well-behaved physical solutions (see also Ref. In the 

case of the plane gravitational wave we have studied in this paper we note that the solution 
is clearly non-perturbative in the coupling constant (3 and still the full solution physically 
makes sense. The "runaway" or "unphysical" solution [Iq does not appear in Eq. ( p.l6| )) is 
naturally discarded here. Besides, we have seen that this kind of gravitational waves will be 



a solution of more general Lagrangian than ( |1 . 1| ) with field equations containing derivatives 
of order higher than the fourth. 

Note also that had we considered the possibility of /5 > in Eq. ( p.l4| ), we would have ob- 
tained a profile function /(p) proportional to the Bessel functions of the first kind, Jo(p/ V^) 
and Yo{p/^/J3). The problem with these solutions is that now the large p behavior is oscil- 
latory with a slowly decreasing amplitude [as {p/^/P)^^^"^-] This destroys the "Newtonian" 
limit of the solution, i. e. the Aichelburg-Sexl profile. 
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note added: After completion of this work we found paper |^ which deals with sourceless 
pp-waves in higher order theories. Our results appear to be completely compatible with 



those of Ref. [p5 |. 



APPENDIX A: SCATTERING OF FIELDS BY SHOCK WAVES 

A further interesting application of the gravitational shock waves was stressed by 't Hooft 
T6| . These geometries are relevant for ultra-high-energy scattering processes (i.e. Planck 
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scale scattering). Let us consider a frame of reference where particle 1 is practically at rest 
and can be described, in a semiclassical approach, by a scalar field satisfying the Klein- 
Gordon equation in the curved background generated by particle 2, which carries Planckian 
energies in our chosen system of reference. In this case one can compute the form of the 
scattering matrix and prove that there will not be particle production of the scalar (or 



other spin) fields. For the axial symmetric case the S-matrix takes the following form 



5(^1,1,^1,2, ^^) 



1 

27r Jo 



M\ P±,i -P±,2 I p) e 



i<^/(p)/4 



p dp 



(Al) 



In Ref. |^2[ the S-matrix have been computed and studied for several sources. It is evident 
from the form of the profile function ( |2.16| ), that it is difficult to find an exact analytic 
expression for the corresponding S-matrix. To illustrate the use of the above expression we 



shall consider a recent interesting result where it was obtained the profile function for 
the ultrarelativistic Kerr geometry 



/(p) = -8Gpe{p - a) Inp + 8Gpe{a - p) 



— In a + 



a? — 



(A2) 



for the boost along the axis of symmetry and where G is the step function. It is still difficult 
to perform the integration in ( [Al| ). We thus approximate the profile function to be constant 
inside the "ring" of radius a and leave the exact logarithmic dependence outside it, as shown 
in Fig. 1^. One can easily check that this approximation is good up to order for small 
p while is also good near p = a where df /dp diverges. We thus decompose the integration 
interval into two pieces for the dimensionless variable p = p/a, form to 1 and from 1 to 
infinity. The final result is given by ISF 



\2~is 



S{q, 



27iq 



2-*^e*^/Vi(g) + q'' (zs Jo(g)5_,,,_i(g) - Ji(g)5i_,,,o(g)) 



(A3) 



where q = {\ pl,i — Rl,2 |)/a, s = 2Gpuj, the Mandelstam variable, and S^^y{z) is a Lommel 
function. A relevant point to see in the structure of the S-matrix is whether it has poles that 



would eventually correspond to bound states of the system. One can see that the above 
expression has no poles in the s variable. This has to do with the fact the source has an 
extended nature, as opposed to the boosted Schwarzschild geometry (i.e. Aichelburg-Sexl 
metric), for which there are poles at is a natural number. 

The same conclusions apply to the wave ( |2.16| ). In fact, now the corrections due to the 
higher derivative theory produces a profile function (that is what matters to the scattering of 
scalar fields) which is finite at p = 0, as we see in Fig. |1|. The same approximation procedure 
can be also applied to the profile ( 2.16| ). This time we approximate the profile for, let us 
say, a; < 3 by a straight line given by f{x) = 0.43a; + 0.116, and the match to the general 
relativistic logarithmic behavior for a; > 3. The non-appearance of poles in the S-matrix, 
again can be traced back to the fact that the /9-dependent term can be seen as an extended 
source (see Eq. (|3.5|) ). 
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FIGURES 



FIG. 1. We compare here the profile function for the general relativistic solution, given by the 
logarithmic behavior, with the profile solution of the quadratic theory. It is evident the regularizing 
effect of the quadratic terms near the particle, located at (p = 0). We also observe that after a 
few mass distances (m = 1/-^/— ^) from the origin the two curves become undistingishables. Here 
X = pI yj—fi, and we have substracted an irrelevant constant from Eq. (|2.16|) , and set to one the 
factor —8Gp. We thus plot Kq{x) + \n{x). 

FIG. 2. The profile function /(p) = f{p) + 8Gpln(a/ po) for the ultrarelativistic Kerr geometry 
and our approximation for p < a to compute the scattering matrix, S. For p > a we take the exact 
(In p) behavior. The approximation is very good near p = and p = a, and allow us to get the 
relevant features of S. 
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